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A GENERALIZED VORTEX THEORY OP THE SCREW PROPELLER 

AND ITS APPLICATION 
By Hans. Reissner 

SUMMARY 

The vortex theory as presented "by the author in ear- 
lier papers has been extended, to permit the solution of 
the following problems t 

(1) The investigation of the relation between thrust 

and torque distribution and energy loss as 
given -by the induction of helical vortex- 
sheets and by the parasite drag; 

(2) The checking of tho thdor em of Bet 2 of the rigidly 

behaving helical vortex sheet of minimum in- 
duced energy loss; . .. 

•(3) The extension of the theory of tho scrow propel- 
ler of minimum energy loss for tho inclusion 
of parasit o-drag distribution along the blades, 

A simple system of diagrams has been dovolo ped -to 
systematize the design of airplane propellers for a~wTdd 
range of practical application. Several typical diagrams 
are presented to- illustrate the method. 



INTRODUCTION 



In oarlior papers (references 1 and' 2), the author 
presented a theory that followed tho conception presented 
by Botz and Prandtl in 1919 (reference 3) of the helical 
vortex sheets loft behind by the trailing edges of tho 
propeller blades. y . ■ — - - 

The relation between tho thrust and the torque dis- 
tribution and energy loss have boon given in a somewhat 
now and more precise form by rational superpos it ion "of tho 
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influence of potential flow and of friction. The theory 
of 3otz was proved more rigorously by extending a method 
of proof used "by Fuchs (roforence 4) for the straight- 
moving airfoil to the rotating air-f-oil. 

A new theorem for the screw propeller of induced and 
parasite drag energy loss can he stated by a very simple 
formula giving the induced inflow angle change of miTTimum 
total resistance (or onorgy loss) as affected "by an. arbi- 
trary distribution of parasite drag along the blade. 

A method of computing the induction of rigid or non- 
rigid vortex shoots entirely differont from the method of 
Goldstein has been followed* The results of this method, 
as far as thoy correspond te rigid shoots and no parasite 
drag, shall bo compared with the results of Goldstein (ref- 
erence 5). In order to derive a system of tables and 
curves for the range of engine powers, speeds, air densi- 
ties, and relative inflow for the design of airplane pro- 
pellers, a great amount of work is necessary and all the 
charts have not been complotod. It was necessary t~o — de- 
velop the circulation function upon which all the blade 
dimensions depend into a power sorios along the blados, 
satisfying certain boundary conditions at the tip and tho 
root. Tho values of the constants were determined by mak- 
ing a mininnim the moan square of the deviation from the 
prescribed inflow angle change of minimum total energy loss. 
Such a mean square can bo defined in different ways, whioh 
either emphasize the correctness noar the tip or near the 
root . 

About 200 tabulation sheets and 86 diagrams have boon 
prepared, embracing tip speed-advance ratios from 1 to 5 
and blade numbers from 2 to 4. Theso diagrams allow im«» 
mediate application for tho ratios 3 to 5 but, for the 
ratios 1 and 2 giving high flight velocity, moro accurate 
values of tho Bessel functions appearing in tho intograls 
will have to bo determined for the nriddle part of tho 
blados. (Dwenty-six of the diagrams are reproduced as sam- 
ple sheets in this paper. 

The work was carried out at Armour Institute of Tech- 
nology with tho financial assistance of the National 
Advisory Committee for Aeronautics* 
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torquo acting on propeller either directly or 
through, a gear from the engine. 

thrust of propeller. 

total, loss of energy (Qu) - Tw) . 

radial velocity. 

angular velocity of propeller . 

propeller speed, revolutions per minute. 

axial velocity of air relative to plane of rota- 
tion of propeller. 



resultant inflow velocity ( V ws + asrs ) . 
induced velocity, 
radius of blade element, 
radius of "blade tip. 



density of medium (air) in which propeller is 
working. 

"blade width. 

number of "blades. ■ ' — ' 



dimons ionless radial abscissa or reciprocal 
advance per unit of circumferential travel. 

angle with piano of rotation of undisturbed^ 

relative resultant inflow (eotan -1 p ; p~ = cot a) 

lift coefficient. 



parasite-drag coefficient. 

torque coefficient ( S-^ 

\,|j.w 5 



) 



thrust coefficient ( -^if V 
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C-g , energy-loss c oef f Iclent ^ 
axial coordinate. 



T, circulation. 

r*, circulation coefficient ^ r ^g- ^ • 

8, angle of parasite drag (6 = cotan -1 C L /C D ). 

i, change of— angle a of undisturbed inflow pro- 
duced by the induction of the helical rortex 
sheets. 

FUNDAMENTAL PROPOSITIONS 

A dimens ionles s radial number 

p = o> r/w (1) 
According to figure 1: 

dL = C L m- dA (2) 

2 

V 2 

dD = C D \x dA (3) 

2 

E = Q, UJ - Tw (4) 

Considering (figs. 1 and 2) that the re-sultant ele- 
mentary lift dL acting on the elementary area dA = bdr 
of a blade is perpendicular to the incidence line of. the 
flow as it is changed by the induced angle i and the 
drag angle 8, the thrust, the torque, and the energy 
loss can be expressed by the following integrals: 

q = X / dL sin (a + i + S) r- (5) 



T = t / dL cos (a + i + 8) (6) 

'•• 0 



E = I / dL w Sin (1 * 6) 

J 0 sin a 



(4a) 
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3?or the sections used in propollor design and for tho 
practical anglos of attack "between the direction of no 
lift and the relative inflow direction, the angle of para- 
site drag changes slowly going from tho maximum angles 
used toward tho zero angle of attack "but increases rapidly 
in tho other direction toward tho angle of maximum lift. 

The lift dL ,. ■ that the air current oxorts on a blade 
element cdr shall be expressed by the circulation r 
produced by tho helical vortex sheets originating at tho 
trailing edges of all the blades in the following manner: 
Bquat ion , (-2) bocomes in this case . . 

dL = C L & V 2 cdr (2a) 

where 

V 2 = u)Sr 2 + w s = w 2 (p2 + 1) = w» (1 + cot^a) = (7) 

sin 3 a 

and where, following Vies elsberger 1 s relation, the lift 
coefficient C L shall be split in the ideal (friction- 
less) part and the deduct ion -caused by the parasite drag. 
This simple superposition of the two effects of potential 
flow and of friction has been proved to bo sufficiently 
accurate and very useful by experience and practice. Put 

C T = 0 T (1 - k6) (8) 
i> L,o 

with the experimental constant k having a value of about 
7.5. ~ 

Those rolations must now be compared with tho Kutta- 
Joukowski-Prandt 1 relation between circulation T and lift 

dL 0 = p,r V dr (9) 

where dL 0 is tho "ideal" lift clement; that- is, tho lift 
element for an airfoil without parasite drag.. 

Comparing now equations (2a), (8), and (9), obtain 

6l,o c 7= 2 r 



(i - ks) = n (i - i=6) do) 

w sin a 



and 

dL = jj, r V dr 
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Introducing thon a dinonsionloss circulation coef- 
ficient by putting 

r = r* ~ (11) 

remembering that, fron equation (1), 



sin a = - - : ■ — t cos a, = ^ . — 

j i + p 3 y i + p 3 

and taking int-a account that JP* and i + 6 aro snail 
quantities of which products and power only up to tho 
second power noed be retained, it is possible to write the 
relations (5), (6), and (4a) in dinonsionloss forn, as 
follows i 

Dimensionless force and energy coefficients may- first 
be introduced: 

H 



Orp 



Qo) 3 

|0, w 

\x\-i 4 

Bo) 2 
jjlw 5 



(12) 

(13) 
(14) 



Then the relations (5), (6), and (4a) appear in the fol- 
lowing pure— number expressions: 

Cq = r dpp r*[l-+ p(i + 8)](1 t W) (5a) 
" Pi 



P 



a 



Orp = j* dp r* Cp - (i + 8)](1 - k6) (6a) 



Pi 
P a 

C 3 = dp r* (i— + 8)(1 + p 3 )(l - k5) (4b) 

Pi 
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In references 1 and 2, three important, problems of 
the Tor t ox- shoot theory of the. propeller are stated: 



To determine the distribution of induced-angle 
changes i of the. inflow along the "blade 
radii produced by given distributions of cir- 
culation coefficient r* (or vortox-dis'tri- 
bution coefficient dr*/dp). 

To determine the distribution of circulation co- 
efficient r* following from a given distri- 
bution of blade width, blade section, and 
blade-angle sotting. 

To determine those distributions of circulation 
coefficient T* and tip radial numbers 
p a = cur a /w that give- the minimum onergy-loss 
coefficient Cj; if the value either of the 
torque coefficient Cq or of the thrust coef- 
ficient Crp is prescribed. . 

In references 1 and 2, methods havo boon indicated to 
solve these problems. *2his paper shall consider only prob- 
lems (a) and (c). It slir.ll bring a new. proof and an exten- 
sion of the Botz theorem of the rigid vortex sheet and shall 
work out one of the indicated methods as described in de- 
tail in the following sections • 



IiHTLOtf ANGrLE CHANG-S IHDUCEIOIT i 20E. MI1TIMTJM E1TSRC-Y LOSS 



"Before -the discussion proper, a ronark shall be 
made about a reciprocal theorem in potential theory and 
its application on the determination of distributions of 
ninimum induced resistance. 

R. Fuchs has proved, by moans of this reciprocal 

theorem, M. Munk ' s theorem of the elliptic distribution of 
circulation and of the rigidly moving vortex shoot for air- 
foils in translation and has suggested that it also might 
bo applied to the problen of airfoils rotating in a relative 
air current and, in this way, the induced resistance x>v 
energy loss of the rigid helical vortex Shoot investigated. 
(See reference 4.) 



namely , 
(a) 

(b) 
' (c) 



The follo\fing discussion shows that tho reciprocal 
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theorem leads to the result that the generalization of 
Munk 1 s theorem to th^ helical vortex sheet is -justified 
up to small quantities of the second order and can even 
"be further generalized to include the parasite drag.' 

The reciprocal theorem (due to Gauss) will first be 
given in the language of potential theory and then trans- 
lated into the language of hydrodynamics. 

Given a closed surface A and two potential func- 
tions Cp^ and cp 2 continuous with continuous derivatives, 
the reciprocal theorem of Gauss can then he pronounced in 
the form 

/dA cp, ^ = / dA cp 3 ^ • (15) 

where the derivatives are to he taken in respect to the 
normal of the surface taken everywhere to the inside (or 
everywhere to the outside). 

If now (fig. 2) this surface, is given "by two of the 
helical vortex sheets, two plane sect±o.ns parallel to the 
plane of rotation and a cylindrical surface coaxial to the 
axis of rotation at a sufficiently great distance from it, 
then only the derivatives at the vortex sheets give a con- 
tribution to the integral "because : By the supposition of 
helica.1 symmetry, the terms of the two sectional planes 
cancel each other; the terms on the remote cylindrical sur- 
face disappear, decreasing more quickly than the recipro- 
cal square of the radial distance; and the terms on the 
helical surfaces also cancel e-ach other because, on account 
of helical symmetry, the potentials are equal and the de- 
rivatives are of opposit-e sign on account of the definition 
of the positive sign of the normal. 

Observing then that 

dA =■ dz dr — 
w 



= dz dr V 1 + p 

and that dqp/dn is the velocity induced by the vortex 
sheet, that is, (see fig. 1) 



|SB * v* = iV = iw J 1 + p 2 
on 
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the reciprocal theorem can be written for this special 
case in the form-. " 

/ r* x i s (1 + p 3 )dp = / r* 3 i x (l + p 2 )dp (15a) 

the integral, taken only over the discontinuities of one 
helical discontinuity , (or vortex) sheet. 

Comparing with (4t>) and assuming for a moment_thc 
paras ite~drag . ailgle 8 = 0, this relation gives g. recip- 
rocal theorem for the energy loss corresponding to one 
kind of circulation function (or T* x written ~dimen- 

sionless) and its corresponding flow deviation i x and 
another kind of circulation Vz (T*s) and its correspond 
ing flow deviation i 3 : namely, ' --- 

G E>1,2 = °E,2,1 ( 15b > 

The idea of. the proof of Fuchs is to consider first the 
energy loss produced by the difference of two circulation- 
deviation distributions T* x , i x » 37*3. i'ai namely, 

c E,l-2 = Ap(F*i ~ I*a)(i x " + P 3 ) 

= °E,1 + C E,2 - /ip(r* x i 3 + + P 3 ) .... 

and, according to (15a) , 

°B,i-a = c e,i + °e,2 - 2 / 4 p^ + p a ^ r** 1 ! 

Nov; for a rigid helical sheet, X being an unknown con- 



stant, the induced velocity V* and the induced "an gTe 
i r would bo. . 



V* 

- = i r =X — ~ - (16) 

1 + p 



Assuming i r . = i^ and inserting this value into the last 

equation of the enorgy loss produced by the differences, 
obtain 

°E,l-2 = C E,1 + G E,2 - 2 */ d PP (17) 

As cither the value of torque Q or of thrust T must 
be proscribed (which one does not nako any difference) , 
prescribe the torque of the engine according, to (5a) and, 
interested for the monont only in the energy loss by the 
induced resistance but not by drag, the condition of keep- 
ing the torque constant is; 
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Pa . Pa. 

= P dpp r* x (i + pi x ) = r dpp r* 3 (i + pi 8 ) 

'i± ■ "pi 

and, if again for ± x tho value appertaining to the rigid 
vortex sheet is introduced, then: 

/dpp r*i(l f A -)= /dpp r*s(l + pis) 

^ 1 + p /■ 

or 

A/d PP r*z = a /dpp r* x + A/d P p s (r*, Xp - r* 2 i 2 ) 

v i + p 3 

Inserting this expression in (17), this equation loads to 



°E 



,2 = g e,i-2 + °h,i + 2A/d PP 3 (r\-i£-- - T* 3 i S ) da) 

x 1 + p 1 * ' 



Nov/ tho last integral is a snail quantity of tho third 
ordor, if T*i i, and A arc all of tho first ordor of 
snail quantities, and therefore 

C E,2 = °E,1 + °E,l-2 < 18a > 

and Betz^ generalization of Munk ' s theorem, to attain the 
minimum induced energy loss for a given torque or a givon 
thrust, is provod to be correct "by the method of 3?uchs. 

3?his significance of (18a) results from the fact that 
tho energy loss 0_ .. ■ must "be positivo and therefore 

Ei , 1 ** & 

G E,1 nrnst~be smaller than any other possible loss C-g^g 
of energy. 

In ordor now to extend the minimum condition to tho 
total-enorgy loss including the -parasite drag given by tho 
drag angle 8, if either torque or thrust of tho propel- 
ler (also including the effect of the drag) is proscribed, 
the energy loss corresponding to the differonco of two 
circulation functions r*, f F*s and to the difference of 
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the corresponding deviation angles i must "be expressed 
in the general form containing the drag angle 6. 

This expression is given by 



Ce,i-2 = /M* + p*)(t\ - r* a ) p x + a - (i 3 + 8)3 (1 - ks) 



= c e,i +c e,2- /ip(i+p 2 ) [r* 1 (i a +0+ r* a ( *i+«)3 d- k6 ) 

Sow the potential part of the flow satisfies again the con 
dition (15a) 

/d P (i + 9 s )r* x i 2 =/d P (i + p s )r* 2 i x 

Therefore, neglecting small quantities of the third order 
like r*iS against the quantities of the second - orcLef 
like r*i' or r*6 : 

C E,1~2 = °E,1 + OB, 3" 2 / d P( 1+ P 3 ) r * 3 (ii + §) - /dp ( 1 + p*)r* x 
Putting then 

l x + -I = X 2— (16a) 

2 ' 1 + p 

the last equation becomes 

G E,l-2 = °E,1 + OB, 3" 2.X/ dp p T* a - /dp(l + P s) r*!6 (18b) 



The third term on the right-hand side can now be trans- 
formed by the condition of prescribed torque, according 
t o (5a) : 

c q = /*pp r*,[i + P (x r -^ y + |)]d- w) 

= /dpp T* 2 '1 + p (i s + 6)3 (1 - k8) 

This equation, if small quantities of higher order than 
the second are again neglectod, furnishes; 



2 x/dpp r* 3 = 2 x/dpp 



On the other hand, 
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20i, J. = 2-7. dp r* x (l 4- p ») ( X 2 + £ ) 

^ 1 + p 2 2 / 

= 2x/d P p r* 1 + 2 /dp r* x (i + p s ) § 

Putting this expression into (l8t>), thore arises finally 
again tho result 

°E,2 = °E,1 + °E,l-2 (l8a) 

which result signifies an extension of Beta's theorem from 
ideal propellers to the actual propeller having parasite 
dra g as follows"; 

The pr opeller--of minimum loss of energy, if torque 
or thrust is -prcscr ibod . must loavo behind its 
trailing blade odgos vortex- she ot s that Will pro - 
duce perpendicular to then induced inf low-angle, 
changes (i) equal to tho difference between the 

p 

inflow-angle changes ^ ■ of rigid helical 

1 + p 2 

surfaces and the half_of tho angles 6 of para- 
site drag of- the blade sections * 

In this theorem, tho factor A remains indctorminato 
unt-11 the relation between the circulation T* and tho 
inducod angle deviation i (or induced velocity V* = iV) 
has been devoloped; whereas the drag angle 5, although 
not varying much for the most important range of angles of 
attack, must bo chosen from experiment according to tho 
shape and tho smoothness of oach section of the blade. 



DETERMINATION OF DISTRIBUTION OF CIRCULATION 



(THRUST AND TORQUE) 



In order to apply tho rosults arrived at in the fore- 
going pages, it is nocossary to answer tho question: 
Which circulation distribution r* would be nocessary to 
produce a desired induced change i of inflow angle? 

This question has not bo answered directly until now, 
oxcopt for tho inflow angle change of the rigid vortex 
sheet by S. Goldstein (reference 5), but the inverse prob- 
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lera, to determine i from an arbitrarily given circula- 
tion r*i is in principle solved by the Biot-Savart in- 
tegral. 

In references 1 and. 2, there has "been derived for 
this purpose a formula different from the Biot-Savart 
integral, which avoids the three-dimensional space inte- 
grations and the inconveniences of apparent singularities 
in the Biot-Savart solution and gives the r resultant in- 
duced angle change i .as the sum of three terms. 

The first term, which is independent of the number 
of blades (that is, for an infinite number of infinitely' 
small blades), signifies the mean value of i on each 
arc between each pair of I blades and can be called the 
self-induction i s of the circulation at a blade point 
on this point. 

The second term represents the 1 effect i^ of the 
circulation existing on the radius inside of a point, on 
this point, using an infinite series of integrals contain- 
ing the circulation distribution in this interval and de- 
pendent on the number of blades. 

The third term in an analogous manner represents the 
effect i e on this point of the circulation existing on 

the part of the radius outside of it. - 



Introducing the Bessel functions of imaginary argu- 
ment of the first and second kinds I n ^(ntp) and K n ^(ntp) 

written in Watson's notation as real functions and writing 
for abbreviation 

nl = m, Ini (nip) a I m V K n \(nlp) « K m 

the formula may be written: 

f* 



i(p) = 



4tt p 2np 



nP • .(20) 



This formula covers the cas.e where the induced velocities 
V* or the angles i remain finite at the axis and at 
infinite radial distance. 



Anyhow it may be useful for applications to give also 
the formula for the case where a hub body, as in figure 
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(2a) makes the radial velocity u zero f-or , r = r^ , that 
is, for p = pi. 

In equations (4), (?), and (9) of reference 1, the 
two integration constants f n ^ 1 ^ and fn^ 2 ^ have to be 
determined for the assumption of a huh "body, so that: 

(a) for p = oo, it is required that u (°) + uC 1 ) = 0, 

+ v ( x ) = 0, W C°) + wC 1 ) = 0, and this 
relation is again accomplished by putting 



'J o. 



(b) for p =r Pi , it is required that u(°) +-u^ 1 ^= 0, 
v (o) + v (i) finite, w(°) + finite and 

this expression brings with it 



f n ( c >= o, (*H) = o, (r*) = o 

V dp / p=Pi P=Pi 
These conditions change formula (20) into 



(21) 



i(p) - 



p* 



4tt p 2n p 



+ Im(p) 



(20a) 



provided that conditions (21) are observed for the choice 
of the function F*. 



This fact would imply that, for every change of~ hub 
diameter (change of p^} , another function V* would have 

to be employed. Yet, if p^ is sufficiently small in re- 
lation to p a and if only T* and dr*/dp are suffi- 
ciently email at the point , then conditions (21) can 
be neglected in view of the fact that, on account of the 
fillets between the hub spinner and the blade, conditions 
(21) are actually only approximately valid. 

If now the problem is advanced to design the I 
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"blades with such widths, cross sections, and angle posi- 
tions that-, for a given speed of the airplane w (cor- 
rected for the influence of hull), an angular velocity of 
the shaft .<ju (speed of the engine), and a given torque 
Q,, thrust, is a maximum or, what is the same, "the energy" 
loss is a minimum, then the distribution of circulation 
r* must "be chosen according to the law (16a) for the in- 
flow angle deviation i previously derived: namely, 




(16a) 



where the 5 l s are the drag angles of the cross sections 
of the Dlade, generally increasing with thickness, that is , 
from the tip to the root of the hlado; and X is a con- 
stant to he deterniinod to give the roquirod torque Q ac- 
cording to (5a) * 



The r elat ion . ( 16a) must then "be introduced into (20a) 
so that the following determining equation for T* arises: 



Pi 

Pa 

p 

This equation looks like a nonhomogeneous integral equation 
of the second kind for the special case where, under the 
intogral, the function to "bo determined appear s * as ' a "deriv- 
ative. Trying to differentiate this equation in such a 
manner that the left-hand sido is transformed into - " 
d / a r* \ 

- — f p— — ) and then introducing this expression as a now 

dp \ r dp J - — - 

function y to arrive at the customary form of the inte- 
gral • equation , one fines that, following a proporty of the 
Bessel function, the series splits into two divergent 
series. It has boon- proved, however , that the screes is 
convergent before this • trans formation , so that this trans- 
formation does not lead to a useful meth.d. It shows only 
that the induced velocity or inf 1 ow-angla change i can 
bo determined as a summation of the vortex " distribution 
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d r*/dp but not as the summation of tho of foots of the cir- 
culation (r*) distribution. 

Tho solution of (22) has boon effected in the follow- 
ing v/ay: Assuming for r* a polynomial series in p with, 
at first, arbitrary coefficients but satisfying tho condi- 
tion that each polynomial torm T*p is zero at the inner 
and tho outor boundaries p^ and p a and that the deriva- 
tive dr*/dp is zero at p = p^ , it is assumed for the 
special case p^^ = 0: 

T* = P 2 Pl Ap r* p , where F * =; p P+s (p a - p) (23) 
p=o 

The angle deviations i produced by this circulation 
function have been computed by (20a), as will be describod 
in detail, in the form 

i = 2 A p i p (23a) 
P=o 

whore the ip' s arG wow known as functions of p and p a . 

Comparing now (16a) and (23a), the constants are to 
be determined by tho equation 

^ A p i p =6v 2 1) (24) 

p=o ^ 1 + p a 2 /p a >p> Pi 

with the condition f-or prescribed t-orquo 




1 + o 3 2/ 



(1 - k6) (5b) 



= / dp P 2 (a p r* p ) |i + , , 3 + 
'i 

This integral is easily calculated as soon as tho function 
6 is prescribed and- gives the value of tho constant \ , 

In order to - find the best representation of tho right 
sido of (24) for a limited number of terms p f tho mean 
square of th© deviation from tho ideal functicn, that is, 
from tho right side of (24) , shall be made a minimum in 
such a way that fro tho doviations from their true values 
is allotted groater weight toward tho tip of the blado, 
increasing with tho area covorod in rotation. 
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This method furnishes the following condition for the 

P=Pi 



coefficients Ap 



dpp X — e_ - f - 2 A p i (p) 



minimum 



Pi 

The differentiations in respect to each of all the coeffi- 
cients Ap gives the following set of equations for the 



A p »s: 



2 A^ / dpp i a (p) ip(p) = X / d P P i n (p) 
P=o * J • / 1 + p* 

dpp S(p) i n (p) (24al 




P 



In each of these equations appears one angle i n . The num- 
ber of these equations is equal to 'the number p x of poly- 
nomials r*p chosen to design the blade "because a runs 
from n = 0 to n = p x . 



The functions i p and i n found by (20a) and (24a) 
express all the constants Ap "by the constant X • The 
constant X must then ho. determined by means of the pre- 
scribed value of the torque coefficient Cq, that is, by 
introducing the Ap , X relations from (24a) into (5b) , 
which then becomes a quadratic equation for X. Of the two 
roots of this equation, only the small and pos'iTive" one has 
a physical meaning. 



CALCULATION OF ANGLES i p INDUCED BY 
THE CIRCULATION TERMS r* p 



Eor the procedure described, it is necessary to com- 
pute. the functions ip(p) by means of (20a). As the 
Bessel functions of imaginary argument I nl (nip) * I ra and 
K n l(nlp) = E m are 'not tabulated as far as they are needed 

here, -it had already been suggested by the author (refer- 
ences 1 and 2) that the asymptotic expros sions~~~bf Nicholson 
(reference 6) be used. Eor the purposo of this " "pr obi em, 
the first terms of these expressions are sufficiently ac- 
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curate down to a value of a radial distance number p = 
0.6; which, for nl = 2, gives an error of 2 percent and, 
for the next terms .with nt = 4 and nt = 6, gives 
errors less than 0.1 percent. In a ' supplementary section, 
it is shown how values for points nearer the origin- are 
determined "by the classical power series of the Bessel 
functions and how the behavior of- the curves, to he con- 
putod, approaching the origin, is ascertained. 

The Nicholson formulas for the Bessel functions of 
imaginary argument as used in references 1 and 2 have heen 
written as follows: 

Introducing as a now independent variable 

u = J 1 + p», v = J 1 + I s , | = J v 3 - 1 (25) 
then Nicholson's formula can he written: 



:„(„> =_L__L (,*/«ZI) 



m 



~m 



and, introducing the development (23) for T*» 

r* p = p p+s (p a - p) 

- (e P +1 f(P + 2) 3 p a - (p + 3) 3 i } 

the term i p corresponding to the contrr ibut i on of r*p 
can be written (using (25)): 



1 P+ 1 / \ 



m+n 



1*1 

4TT 



.017 (▼-!) 



ai-p 



(p+2) 2 Pa - (p+3)Vv 3 - 1 



rau 



\ u+l ) 



a/ u(Tf -1) 



Pv 



'^Pi 

JL+P_ 
-mv (7-1) a 

EC- p 



(P+ 3) 8 Pa " (P + 3 )V v2 - 1 1} 



(20b) 



*5 
o 



d> 
o 

S3" 

P 

H" 

O 



!2j 

o 

ct- 
(P 

o 



-0 
CJl 
O 
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In expression (20b) , thore appear- four integrals. 
Each, of the integrands has "been computod for tho valuos 
of m = 2., 4, 6, 8, and 10 and for values of p = 0, 1 ? 
and 2 . '• 

The calculation showod, as is to to seon 1 in figures 
3, 4, and 5, that tho integrands increased so rapidly 
with the abscissa v for I m and v a - v for K m that 
it was necessary— to change tho scale very often in plot- 
ting them along tho abscissa. 

This change v/as particularly necessary in order to 
bo r.blc to draw the intogral-. curves , that is, in ordor to 
detormine tho areas for a sufficiently great number of 

abscissa valuos | = ^ v 3 - 1 or v = «/| 2 + 1. 

Each of tho integrals in (20b) was represented by 
two curves corro spending- to the tiro terms under the inte- 
gral signs in (20b) from whicn tho different curves re- 
sulted. 

In this way for the induction from the inside, the 
sets of curves for the different tip radii numbers 
u a = 2, 3, 4, and 5 in figures 6 and 7 and, for the in- 
duction from the outside, the sets of curves in figures 
8 and 9 were drawn for the different tip radii numbers 
u a , the different orders in p, and the different orders 
in m. 



For the induction at each one of the points u =vp 2 + 1 
between u = 1 and u ■= u a , the integral as well as 
from the outside as from tho inside had to be multiplied 
by tho value of K m (p) and Im(p)i respectively, at the 
point considered, as equation (20b) shows. 

This task was rather laborious by tabulation as well 
as by drawing, but nothing simpler could bo found. 

Having then determined tho angle inflow induction ip 
from inside and outside for a certain tip radius number 

u a = J p a 2 + 1 and a number of points between u = 1 (p = 0) 

and u a , both inductions were superposed in figures 10 to 
16 so that a graphical picture of the total induction for 
every p and for every m was gained. 

Finally, for oach p separately but for all tho 
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values of m that have been computed, the induction curves 
have "been a<ided and a remainder of the m series has "been 
added "by utilizing the fact proved elsewhere that the con- 
vergence of the aeries in m approaches the series 

2— with sufficient approximation. 



These resulting curves are shown in figuros 1? to 24. 
All these curves do not contain the first term of "self" 
induction on the right side of (20b), which is, vory 
simply, to he calculated and added to the "distant 11 induc- 
tion curves. 

It may ho remarked here that, for reasons of "being 
able to compute (20b) graphically, the functions chosen 
in (23) have no vertical tangent at the point p = pa» as 
is the caso for elliptical distribution. 



finite. Tot the inclination of the tangent to the curve 
r* can bo made as groat as desired by ovon a singlo term 
of (23) if the condition of minimum energy loss would de- 
mand it. . . ._ 

]?or the following computation,- two blades (l = 2) 
have been assumed but, in the future c omput at ions , many 
of the integrals can bo used for other blade nunfbers as 
only the order number m = nl (n = 1, 2, 3 . . «) is de- 
cisive. 

In references 1 and 2, it has boon shown that tho 
convergence of tho scries in (20b) is of the charactor 
m~" s but it remains to be investigated how the right side 
of (20b) bohavosfor very snail values of mp because "of 
the infinite valuo of tho Bessel functions E of the sec- 
ond kind for mp = 0 and because tho approximate repre- 
sentations of Uicholson for not too small values of the 
argument mp = nlp^ have been used. 

3?or small values of mp, the two Bessel functions 
of imaginary argument can bo written 



m 



Such a proposition would havo made 





m 
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The terms used for the power aeries for r* = Z A p r*p 
are : 

• r* p = p p+3 ( Pa - P ) (23) 

giving for the terms under the integrals in (20) 



• The first integral in (20) , if it is to he computed 
only for a very small upper limit p, is then written: 



Evaluating this expression and multiplying by K m /2Trp, the 
result is: 



Km 
2TTp 



s dr* \ pP+i 

\ d| / 4 tt m 



(P + 2) 



Pa 



p+2+m 



- Cp+ 3) 



p+3+m 



(a) 



This equation proves that the curve showing the induction 
of inflow-angle change i by the vortices inside of the 
point considered starts from zero at p = 0 with the an- 
gle tangent 



tan cp 0 



■ — r 

4ttiu [_ 



PaP 



P (P+D (P+2) 3 _ p+i (p+2) (p+3) 



2 -t 



p+2+m 



p+3+m 



Hence, for all values of p > 0, the induction curve .is 
tangent to the abscissa axis and, for p = 0, 



(.tan qp 0 )__- = 



Pa 



P =0 Trm(m + 2) 

The second integral of (20) must also be considered 
in order to judge the trend of the contribution to the 
angle function i. It may be repeated here 



•J n 



Now the integral has been computed by means of the asymp- 
totic approximations of Nicholson for all lower limits 
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down to = 0.6. Therefore only the additional contri- 

butions from p- to p x are considered: namely, 




Excluding at first the particular cases p+2-m=0 or 
p + 3. - m = 0 and considering that the lowest value of m 
is 2, one sees, that, in the limit p — > 0, the contri- 
bution is zero and therefore, .as the total integral from 
P *° Pa i -. s finite and has to "be multiplied hy the zero 
value of I m /2TTp, the induction curves for the influence 
of the outside vorticos go to zero at p = 0 tangentially 
to the abscissa axis, for m > 2, except the particular 
cases. ' 

The particular oases are: 



for 


P ■ = 


0; 


m 


= 2, 


m = 3 


for 


P = 


l; 


m 


= 4 




for 


P = 


2; 


m 


= 4, 


m = 5 



Galling the exponent that "becomes zero q., the ex- 
pression to be evaluated is: 

a q. 

( Pl q" P ' ) ' o = Pi q lo e Pi - p q lo S P = log 

Now for the orders p = 0 and p = 2, one has to apply 
the preceding result to the first term inside the brackets 
of (b) ; whereas, for p = 1, the second term gives the 
particular case. . 
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Thorc results 



p = 0, 



n 



= 2 




P 

p = 2, n = 4 





16p a log £i - 25( Pl - p) (^ 2>4 ) 



p = 1, 



n 



- 4 




Those oppressions show that, for the lower Unit p = 0, 
the contributions arc zero and therefore, as the contribu- 
tion produced by the vorticos at points with larger values 
of p computed with Nicholson^ formulas are finite in the 
values proper of the integral and zero after multiplication 
by I n /2Trp, it appears thnt the induction of the outside 

vortices on the root point in the particular cases previous- 
ly computed (for p = 0) is also zero. 

Yet the influenco of the logarithmic tern appearing 
in ("^0,2^* ^2,4^* ai1 ^ ^1 4^ causes the tangent of tho 
induction curve to be perpendicular to the abscissa axis 
at p = 0. 

DB TEEM I HA 2 I ON OP THE C OH STAKE A AND OE THE CIRCULATION 



The torque coefficient has been expressed by the inte- 
gral (5b). This integral shows that, to each series term 
of 



FACTORS EROM THE TORQUE COEEEICIENT AND EROM THE 



INDUCTION OE MINIMUM LOSS 



= s A p r* 



r* = p p +2 ( Pa - P ) 



(23) 
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there corresponds a contribution of the torque coefficient 
of minimum loss of energy: ' ~ ~ ■ 



.Pa 



*PP P+3 (pa~P> 



1 + A 



TTp + a p 



(l-kfi) (5c) 



The integration leads to (6 is assumed to "be constant and 
in one place 6/2 is neglected against unity) J 



1-kS 



p+ l 



+ Pa 



P+3 



p+1 ! " (p+2).(p+3) 
+ *Pa / 



Pa 



p+5 



(p+4)(p+5) 



] 



P^ A A P + i 

A PP + A 



+ Pa 



P+5 



l+"p* 
1 



(p+4) (p+5) 



Pa 



Pa 



p+ t 



l+p ; 



2 (p+5) (p+6) 



(5d) 



As the terms for T* have beon used only for the num- 
bers » p = 0, p = 1, and p = 2, there are also only 
throe contributory terms for Cq: namely, 



C Q,0 
l-k8 



Pa + 



Pa 1 



Pa a 
20 



arc tan p a - 



log (1 + p a a ) 



_ 5 6 

Pa , e Pa 



20 



+ 8 



60 



Mb- 



l-k8 



Pa 4 Pa 6 
12 



+ i^ 
SO 



s 



~5"+pa arcten p a - \ log (1 + p a 2 ) J 



+ 5 ^ 
84 



^2 = M 



1-kS 



r Pa 
6 



20 



Pa Pa 
42 + 2 



log(l + p a 3 ) - p a + arctanp'a 
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Tho total torque coefficient is then to to wrltton as 
follows: 

°q = A 0 °Q,0 + A l C Q,1 + A 2 C Q,2 < 5o > 



As not only the C' c 



terms "but also the factors A , 



, QjP _ — — p> 

according to (24a) aro linear functions of X , there will 
finally arise a quadratic equation for 7\. 



Tho oquations determining tho factors Ap have "boon 
stated, in (24a). They "may he treated now sonowhat more 
explicitly in such a way as to show how fcise curves computed 
to give the induction of inflow angle ip by tho corre- 
sponding circulation term T^p are to ho used to find the 
coefficients of the equations (24a) dctornining tho factors 

An. 



Writing then equation (24a) in tho form: 

Ao foo + -^l f oi + A 2 fos = A f 0 - ~ go 

Ao f io + Ai fn + Ag f 13 = \ f 2 - i gi 

^O f©0 + -^-i fsi + &3 ^23 = ^ ^3 ~ — Sz 



(24h) 



the factors f and g have the following significance: 

nPa 



•no 



£>n 



which shows that 



i n ^-p 



fn = 



a pp i 



1+p' 



n 



dpp i n 



(27) 



J 



f np - fpn 
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NUMERICAL RELATION BETffEEII THE FACTOR X OF ILEAL 
ANGLE CHANGE FUNCTION i AND THE PRESCRIBED 
TORQUE (POWER) COEFFICIENT 

As a first case has "been troatod the propeller of tip 

radius nunher p a = */ 34 = 4.9. Formula (5d) gives 

C Q,,0 - " k8)(141.0 4- 125. 9X + 230.48) 

0 = (1 - k6)(45l.O + 419. 9X + 806.05) 

C Q,2 = ( x ~ k8)(l, 613.0 + 1,489. 7 X + 2,982.08) 

Supposing now that equations (24h) for A Q > A x , A 3 have 
"been solved in the for:.ii 



A 0 = a 0 X + S 0 S . ._ ■ 

A^ = a J.X + 8 x $ 

A 3 = a 3 X + 0 3 S 

and that the value of the torque (or power) coefficient 
Cq is proscribed hy the perfortaan.ee of the onglno and the 
airplane, then the following quadratic equation for X is 
oht aincd S 

Cq = (1 - k8) [(a 0 X +• 3 0 8)(141.0 + 125. 9X + 230.45) 
+ (a x X + 8 1 S)(461.0 + 419. 9X + 806.08) 
+ (a 3 X + 3 3 8)(1,612.0 + 1,489. 7X + 2, 982.08)] 

or 

X 2 (l25.9ao + 419. 9a x + l,489.7a 3 ) + A[(l41.0+ 230.48)or 0 
+ (461.0 + 806.06)a 1 + (1,612.0 + 2,982.0S)a 3 + 125.988 
+ 419.960! + 1,489. 78B 3 ] + 8 o 8(l41.0 + 230.48) 
+ 8 1 8(461.0 +.806.06) + B 3 5(l, 613.0 + 3,983.06) 
= 0^/(1 ~ k8) 
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Putting for a practical example 6 = 0.02, this equa- 
tion takes the form 

A 3 (125.9a 0 + 419.90,3..+ 1 , 489 . 7a 2 ) • + X(145.6a 0 + 477.1a! 

+ l,671.6a 2 + 2.52£ 0 + 8.463, + 29.793 2 ) + p 0 2.91 + Bi9»51 

+ ^ 0 3S.43 - ^ — = 0 (26) 

2 1 - k8 

NUMERICAL DETERMINATION OE .THE FACTORS A p OE 
CIRCULATION r* AND ANGLE CHANGE i^ 
Case of p a = 4.9 

In equation (24b), the factors fnpp a « ^np a » an ^ 

Sn»p a aro definite intograls that must be determined from 

the angle-change curvos i n _ , which had "been computed 

as the sums of the self-induction, the outsido and the in- 
side induction of the c irculations r*n n each -Doint 

of the radial abscissa p and which are shown in figures 
17 to 24» The integrands of the integrals f and g of 
formulas (27) have boon computed and plotted as figures 25, 
26, and 27. From those curvos, the areas havo been taken. 

Taking these values, the equations that havo to deter- 
mine the factors Ap aro, for the radius number p a = 4»9, 
numerically given by 

A 0 0.1067 + A 1 0.347 + A s 1.236 e — 0.3 - - 0.2153 

P a 

A 0 0.347 + A x 1.331 + A 2 5.269 = 0.859 - - 0.696 

Pa 2 

A 0 1.236 + A - 5.269 + A 3 22.42 = 2.977 - 2.457 

1 Pa m 2 

The solution, which had to be computed very exactly because 
it is rather sensitive to small changes in the coefficients, 
was found to be: 
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A 0 = 



9.4556 — - 1.72385 -| 
Pa 2 



Aj. = -4.0469 



A - 0.23539 -| 



Pa 



A s = 0.56259 — + 0.040708 A 
Pa 2 



(24c) 



so that 



a 



= 1.930 
= -0.8250 



a 3 = 0.115 



B 0 = -0.8619 
B 3. = -0.1177 
B 3 = 0.020354 



The quadratic equation (26) dotormining tho factor X from 
the torque coefficient Cq is then to. he written as fol- 
lows: : 



A 3 68.0 + A 76. 6 - 2.95 - 



1 - k8 



= 0 



(26a) 



To choose a value for tho ( d im ens i onles s ) torque (or power) 
coefficient 0q appropriate for the tip radius number 
(circumferential path per unit of advance) p'q, = 4«9, one 
has to remember that so high a value corresponds to rather 
low power, low speed, and high revolution values', for in- 
stance 



p = 


160 hp. = 


160. 


75 


kg 


-m/soc 


n = 


1,190 r.p. 


n. , 


ui = 


niT 
30 


= 200 sec.*" 1 


h = 
w = 


18,000 ft. 
123.5 in.p. 


h. = 


1 
2 
55 


Uo 

r. 


= -i- kg-m-sec 
16 

/ sec 


d = 


8*54 ft. 











m 160 x 75 x 200* m l5fZ 
H 1/16 x 55 5 

and, with k = 7.5 ind 6 = 0.02, 



1 ~ k6 



0.85 



- 18.0 
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Equation (26a) "becomes 

+ x 0.113 - 0.300 = 0 

As only the positive root has a physical moaning, one ob- 
tains 

X = 0.2225 

It p.ust ho remembered that this value of 7\ is valid if 
the ideally host distribution of induced angle changes is 
obtained. If the best distribution is only approximately 
obtained, the value, of X may chango somewhat from the 
preceding value. 

The idoal distribution of inflow anglo chango was 
stated to bo 

P 6 
iid = X - ~ 

1 + p s 2 

It must now bo shown how closely 'the ideal inducod anglo 
chango function i^ has been approximated with only the 
throe terms corresponding to p = 0 , 1, and 2 by the meth- 
od of the loast avorago absolute .error of equations (24a), 
(24b) , and (27) . 

In the first column of table I are given the computed 
ideal induced angles and, in tho noxt column, the- angles i 
induced by the circulation r* 

r* = a 0 r* 0 + A 1 r* x + a s r* 2 

taken from (23) and multiplied by tho factors A 0 , A x , and 
A s , respectively. 

Tho factors A 0 , A x , and A s can now bo computed 
explicitly after the values of X and" 5 are given from 
( 24c ) , as follows : 

A 0 = 1.93 X 0.2225 - 0.0172 = 0.413 

A x = -0.825 X 0.2225 - 0.0023 = -0.1858 

A 3 = 0.115 x 0.2225 - 0.00041 = 0.0259 
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Table I 



p/Pa 


A 0 io 




A 3 is 


2 A p i p 


s A p r* p 


0.1 


- 


0.109 


-0.02412 


- 


0.350 


.15 


-0.036 


- 


- 


- 


- 


.2 


.01509 


.1042 


-.04605 


0.07324 


.964 


.3 


. 1102 


.0382 


-.0427 


. 1057 


1.970 


.4 


.1818 


-.082 


-.0184 


.0814 


1.85 


.5 


. 2265 


-.1974 


.0373 


.0664 


1.67 


.6 


. 2490 


-.302 


.109 


.0560 


1.56 


.7 


. 2423 


-.379 


.184 


.047 3 " 


1.41 


.8 


.2135 


-.418 


.2466 


.0421 


1.28 


.9 


.154 


-.366 


.248 


.036 


.92 


1.0 


.053 


-.141 


.118 


.031 

■ 


0 



Fijuro 28 shows tho two curvos of inflow angle change 
plotted to the sane scale. The ideal curve for nlnlnuri" 
energy loss agrees quite satisfactorily with tho curve ob- 
tained fron the integration over the throe circulation 
(cr vortex distribution) torus r* 0 , P\ , and V w z down" 
to the third of tho radius; farther inside toward the huh, 
tho curve found fron the circulation power series, al- 
though it shows the sane "behavior in rospoct to a peak 
and tho sudden decrease, needs sone shifting over to tho 
center. This change could he effected by a fourth tern 
of the power serios in r* but, as the higher ^orns shift 
the circulation and also the inflow angle change noro to 
tho tip, it night bo nore satisfactory to add the lower 
tern r* , which had been avoided beforo "~~ T " - 



r*_! = p(pa - p) 

that is, tho sinple parabola. This tern was not taken in 
order to prevent singularity at the center but, as an in- 
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ncr circle of the hub has to ho cut out anyhow, as has Icon 
explained in an earlier section, this singularity would 
not appear and would have no effect. This result proves 
that, with very few terms of the power "series for the cir- 
culation, that is, for the thrust distribution, the ideal 
inflow deviation of minimum loss of energy can bo approxi- 
mated very closely. 

In table I and figure 27 arc given and plottod the 
ordinates for the resultant circulation. 

r* = a 0 ^r* 0 + a x r* x + A a r* 

The curve shows a peculiar sudden docroaso toward the root, 
To the ideal inflow angle chango 

iid =x— e— - £ 

corresponds an (idoal) circulation function r*id> which 
can be approximated the more closely tho more terms r*_ 
are usod. Still a swf f iciontly oloso approximation will 
always bo possible with a small numbor of appropriate sor« 
ios terms. * 

The othor curves cor responding to figure 28 and to 
the othor tip radii numbers p 0 = //3, */H , and </ 15 will 
bo determined in the samo way. 



Armour Institute of Technology, 

Chicago, 111., October 17, 1939. 
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Figs. 7,8 
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Figs. 15,16 
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